Primes Appearing in Tower Factorization of Integers 

Patrick Devlin, Edinah Gnang 
April 23, 2012 



Introduction 

Definition of Tower Factorization 

For each positive integer n, recursively define the tower factorization of n as follows: 
[0] If n = 1, then its tower factorization is just 1. 

[1] If n > 1, let n — p^^ ■ ■ ■p'j^'' be its usual prime factorization. Then the tower factorization 
of n is given by 

where /i is the tower factorization of e^. 



In other words, to obtain the "tower factorization" of n, we first write n into its prime fac- 
torization n = ' ' ' P'k' ■ Then we factor each exponent > 1 into its prime factorization 

Bi —Pii ■ ■ ■ Pi ■ We continue to factor each subsequent exponent of those prime factorizations 
until n is writen in the form 

n = [[P^ 

i 

for some prime numbers pi,pij, .... Then by the fundamental theorem of arithmetic, this tower 
factorization is unique up to reordering. 



Number, n 


Tower Factorization 


Primes q with n G M (q) 


1 


1 





144 


2(2^) . 32 


{2,3} 


725 


5(2^) 


{2,5} 


15379 


7-13^ 


{3,7,13} 


37349 • 11669921875 


133.17.1l(7^-5<3^') 


{2,3,5,7,11,13,17} 



Table 1: Examples of tower factorizations and to which M{q) each number belongs 



Statement of Problem 

Let q be any prime number (e.g., q — 2), and let M{q) be the set of all positive integer^ whose 
tower factorization (as above) contains the prime q at least once. 

Then we wish to determine the density, d{q), of the set M{q). That is, what is 

,^ |M(.)n{i,2,...,iV}| 

N^oo N 

and does such a density even exist? (Stated more colorfully, d{q) is the probability that a positive 
integer 'chosen at random' contains the prime q in its tower factorization.) 

^For later notational convenience, for all q, we set 1 ^ M{q). 
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Exact Formulas for d{q) 



Theorem 1. Let q be any arbitrary prime, and let M{q) be as previoulsy defined. Then the 
density d{q) exists and is given by 



d{q) = 1-1 



1 



n 



1-1 



E - 



P J — P 

^ ' m.eM(q) ^ 



P 7^ <?; prime 

Moreover, if M^iq) is defined M^iq) := {z G {0, 1, . . .} : z ^ M{q)}, th en we also have 



d{q) = 1-(1-- 

q 



n 



p ^ q. prime 
-1 



P / ^ — ' P" 

^ ^ meM<=(q) ^ 



ymeM^iq) 



n 



p. prime 



P J — P 



Proof. Since M{q) C {1,2,3,...} and ^ M'^{q) C {1,2,3,...} convergent exact formula for 
1 — d{q) can be found in the foUowing way, establishing its existence. Let n G {1, . . . ,N} be 
chosen at random (wc will let N go to infinity) . Then for all primes p and all nonnegative integers 
m, let D(p, m) be the event that exactly divides n (that is divides n, but does not). 

Then the probability that n is in M(q) is given by: 

P(?T. e M(q)) = P [{q divides n) U (p™ exactly divides n, where p q, and to € M (q))] 
{q divides n) U I [J [J D{p, m) 

\p 7^ q^ prime m£AI(q) 

These events are not disjoint, which makes this probability difficult to deal with. However, by 
considering ¥{n ^ M{q)) = 1 — P(n S M(q)), we have the simplification 



1 - P(n £ M{q)) 



{q divides ?i) U \^ \^ D{p, m) 

p ^ q, prime m^M{q) 



fl y D{p,m) 

p =^ q, prime \ mgA/(g) 



{q does not divide n) n 
As TV goes to infinity, these events are independent. Therefore, letting N oo, we have 
1 - d{q) = P((? does not divide n) • P Q U 

p ^ prime \ mGM(g') 



^{q does not divide n) ■ 

p ^ q, prime 



1 



y D{p,m) 

KnieM{q} 



Now for different values of m, the events D{p, m) are clearly disjoint. Moreover, as N goes to 
infinity, we have 



{D{p,m)) = P(p™ divides n) ■ P(p'"+ docs not divide n given does) = — '1 



1 



1 



^Note that {0,1} C A/'=(g) 
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Thus, this enables us to simplify our formula to 



d{q) = 1 — W'{q does not divide n) ■ 

P 7^ <?> prime 



= 1-1 



n 



p ^ q, prime 

On the other hand, since we have 



1- 1 



P / — , . P 

' meM{q) ^ 



meM(< 

-) E - 



1 - 



1-Pl (J i:)(p,TO) =P (J D{p,m) 

\0<meM{q) J \0<m^M{q) 

then if we define M^q) = e {0, 1, . . .} : ^; ^ M{q)} (with {0, 1} C M^(g)), we have 



diq) 



1 — P(g' does not divide n) • 

P ¥^ Qi prime 



U ^(f^'H 

^0<m^M(g) 



1 — f[q does not divide n) ■ JJ 

P 7^ *?5 prime 



1- 1 



n 



p 7^ q, prime 
-1 



1 - 



E - 

^6M<!(g) 

E - 



1 - 



P J — , P 



1-1 E 

I meM<=(g) 



n 



1 



-) E - 



' meM'={q) ^ 



p prime 

as desired. 

Proposition 2. If I(n) is the indicator function for the event n e M'^{q), then we have 



i-d{q)= n 



and it is also true that 



p prime 



1 — d{q) = lim 

S->1 + 



p / — ^ p 

^ ^ m=0 ^ 



(i-^)E 



□ 



Proof. This follows because I{n) is multiplicative in that if gcd(a, 6) = 1, then I{ab) = I{a)I{b). 
Then the last theorem can be rewritten as 



1 - d{q) 



-^)^ n 



p =^ q, prime 



n 



p =^ q, prime _ 



-) E - 

1\ g /(pm) 



P / — P 

' ' m=0 ^ 



= n 



P / "'^ — ^ P 



p prime 

Then using that 7(n) is multiplicative, we obtain the second the second formula similarly. □ 
Using these formulas, we can obtain the following bounds. 
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Theorem 3. For each prime p ^ q, let S{p), T{p), A{p), and B{p) he integer subsets such that 
Sip) C M{q) C T{p) and A{p) C Af=(q) C B(j)). Then we have 



dig) > l-i^l 
d{q) < 1 - (^1 
d{g) < 



dig) > 1-1 



n 

p ^ q, prime 

n 

p ^ q, prime 

n 

p ^ q^ prime 

n 

p ^ q, prime 



1-1 



1-1 



- E - 

l\y I 

' teT(p) ^ 

E - 
F - 



beB{p) 



Proposition 4. Le< P &e any set of primes with g ^ P, let S C M(g), let A C M'^iq) C _B. 
r/ien we /ia?;8j 



dig) > i-(i-i)-n 
> i-fi-^)-n 



i-(i--)e- 



< 1 



V pep 
g« - 1 



i--)E^ 

^ ^ beB ^ 



n 



1 



1 Z-^ -n° 



1 aeA ' 

pi 



Proposition 5. For all q, we have 



-7?+i-l Ciq + l) 



gi - 1 Ciq) ■ 



Moreover, we have 



1 1 

< dig) — < 



q 4 • 2"? - 4 • 29 

Corollary 6. The sequence diq) is strictly decreasing, and it decreases to 0. 



^This gives rise to algorithms for computing d(q) 
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Rigorous Computer Bounds for Specific Values of d{q) 



The conjecture d{2) = l/VS is false. 



Using the previous bounds on d{q), the following values were found. The floating point values 
in the fifth column are a rigorous lower bound for d{q) (on top) and a rigorous upper bound for 
d{q) (on bottom). Full values of these bounds are in the file, but they are truncated to just 
35 digits here. 

The columns p, a, and s are parameters for the algorithms used. These parameters correspond 
to the size of P, (roughly) the size of A, and (roughly) the size of S as in Proposition |4l 
Notably, p is the number of primes used in the estimation, which seems to matter much more 
than the size of A [affecting the upper bound] or S [affecting the lower bound] . 



Data and code are available on request. 



q 


P 


a 


s 


Bounds (First 35 Digits [More in TgX File]) 


Digits Known 


2 


25000 


20 


20 


0.577 350 376 056 807 813 001 171 222 749 099 03 
0.577 350 485 047 678 584 952 747 233 500 637 55 


6 


3 


6000 


100 


100 


0.388 807 379 263 994 405 608 000 000 000 000 00 
0.388 807 379 271 511 226 974 000 000 000 000 00 


10 


5 


5000 


100 


100 


0.215 118 984 695 585 620 327 888 615 736 044 88 
0.215 118 984 695 585 620 331 080 414 300 988 94 


19 


7 


2500 


100 


100 


0.146 500 891 228 438 042 819 116 915 103 810 81 
0.146 500 891 228 438 042 819 116 915 105 137 03 


29 


11 


2000 


200 


200 


0.091 134 581 055 674 121 650 272 316 314 808 81 
0.091 134 581 055 674 121 650 272 316 314 808 81 


44 


13 


2000 


200 


200 


0.076 979 810 520 294 777 519 659 200 891 501 69 
0.076 979 810 520 294 777 519 659 200 891 501 69 


52 


17 


2000 


200 


200 


0.058 827 124 602 119 403 676 736 708 884 910 96 
0.058 827 124 602 119 403 676 736 708 884 910 96 


60 


19 


2000 


200 


200 


0.052 632 482 973 467 517 964 355 534 025 063 33 
0.052 632 482 973 467 517 964 355 534 025 063 33 


60 


23 


2000 


200 


200 


0.043 478 317 889 484 083 344 293 667 695 938 90 
0.043 478 317 889 484 083 344 293 667 695 938 90 


60 


29 


2000 


200 


200 


0.034 482 759 519 907 038 838 884 404 979 223 96 
0.034 482 759 519 907 038 838 884 404 979 223 96 


60 


31 


2000 


200 


200 


0.032 258 064 741 450 054 595 016 325 700 965 71 
0.032 258 064 741 450 054 595 016 325 700 965 71 


60 


37 


2000 


200 


200 


0.027 027 027 030 566 683 523 134 092 350 343 62 
0.027 027 027 030 566 683 523 134 092 350 343 62 


60 


41 


2000 


200 


200 


0.052 632 482 973 467 517 964 355 534 025 063 33 
0.052 632 482 973 467 517 964 355 534 025 063 33 


60 


43 


2000 


200 


200 


0.024 390 243 902 660 852 384 118 730 197 418 93 
0.024 390 243 902 660 852 384 118 730 197 418 93 


60 


47 


2000 


200 


200 


0.021 276 595 744 684 328 187 880 387 051 387 64 
0.021 276 595 744 684 328 187 880 387 051 387 64 


60 
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q 


P 


a 


s 


Bounds (First 35 Digits [More in TgX File]) 


Digits Known 


53 


2000 


200 


200 


0.018 867 924 528 301 941 256 223 883 281 009 21 
0.018 867 924 528 301 941 256 223 883 281 009 21 


60 


59 


2000 


200 


200 


0.016 949 152 542 372 882 208 592 895 018 030 94 
0.016 949 152 542 372 882 208 592 895 018 030 94 


60 


61 


2000 


200 


200 


0.016 393 442 622 950 819 885 416 820 821 976 79 
0.016 393 442 622 950 819 885 416 820 821 976 79 


60 


67 


2000 


200 


200 


0.014 925 373 134 328 358 212 292 786 538 441 94 
0.014 925 373 134 328 358 212 292 786 538 441 94 


60 


71 


2000 


200 


200 


0.014 084 507 042 253 521 126 969 339 106 717 70 

0.014 084 507 042 253 521 126 969 339 106 717 70 


60 


73 


2000 


200 


200 


0.013 698 630 136 986 301 369 915 228 058 392 40 
0.013 698 630 136 986 301 369 915 228 058 392 40 


60 


79 


2000 


200 


200 


0.012 658 227 848 101 265 822 785 626 836 554 17 
0.012 658 227 848 101 265 822 785 626 836 554 17 


60 


83 


2000 


200 


200 


0.012 048 192 771 084 337 349 397 641 437 357 10 

0.012 048 192 771 084 337 349 397 641 437 357 10 


60 


89 


2000 


200 


200 


0.011 235 955 056 179 775 280 898 877203 211 62 
0.011 235 955 U5G 179 775 280 898 877 203 211 62 


GO 


97 


2000 


200 


200 


0.010 309 278 350 515 463 917 525 773 198 999 20 
0.010 309 278 350 515 463 917 525 773 198 999 20 


60 


101 


2000 


200 


200 


0.009 900 990 099 009 900 990 099 009 901 185 36 
0.009 900 990 099 009 900 990 099 009 901 185 36 


60 



In this table, perhaps the most important thing to note is that it disproves the conjecture 
d(2) = 1/73 that was strongly suggested by experimental data and early attempts to bound 

d{2). These bounds prove 



0.5773504851 > d{2) > 0.5773503760 > 0.577350269189 . . . = l/Vs. 

It is curious to note how close d{2) is to 1/ \/3; these values agree for 6 digits after the decimal, 
but not the seventh. 
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